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Abstract. We use the Planck HEI channel maps to make an all sky map of yu-distortion fluctuations. 
Our /r-type distortion map is dominated by the y-type distortion contamination from the hot gas in the 
low redshift Universe and we can thus only place upper limits on the /r-type distortion fluctuations. 
Eor the amplitude of yu-type distortions on 10' scales we get the limit on root mean square (rms) value 
< 6.4 X 10“^, a limit 14 times stronger than the COBE-EIRAS (95% confidence) limit on the 
mean of (/i) < 90x 10“^. Using our maps we also place strong upper limits on the auto angular power 
spectrum of /i, and the cross angular power spectrum of /i with the CMB temperature anisotropies. 

The strongest observational limits are on the largest scales, {{i + l)/(2n')C^^|f=2-26 < (2-3 ± 
1.0) X 10“^^ and C(i + l)/(2n')C^^|r=2-26 < (2.6 ± 2.6) x 10“^^ K. Our observational limits can be 
used to constrain new physics which can create spatially varying energy release in the early Universe 
between redshifts 5x10^ < z < 2x10^. We specifically apply our observational results to constrain the 
primordial non-Gaussianity of the local type, when the source of /i-distortion is Silk damping, for very 
squeezed configurations with the wavenumber for the short wavelength mode 46 < ks < lO'^ Mpc“' 
and for the long wavelength mode k]^ w 10“^ Mpc“^ Our limits on the primordial non-Gaussianity 
parameters are /nl < 10^, tnl < 1-4 x 10*^ for ks/^L w 5 x 10^ - 10^. We also give a new derivation 
of the evolution of the /i-distortion fluctuations through the y-distortion era and the recombination 
epoch until today resulting in very simple expressions for the cross and auto power spectra in the 
squeezed limit. We also introduce mixing of Bose-Einstein spectra due to Silk damping and y®^-type 
distortions. The //-type distortion map and masks are now publicly available. 


Keywords: cosmic background radiation, cosmology:theory, early universe 


Contents 


1 Introduction 1 

2 //-type distortions in Planck channels 3 

3 A //-type distortion map from Planck HFI data 4 

4 Constraints on the //-distortion fluctuations 6 

4.1 Amplitude of fluctuations 6 

4.2 Auto and cross power spectrum 7 

5 Theory of fluctuations of //-type distortions from Silk damping 11 

5.1 The auto and cross power spectrum of //-distortion fluctuations 15 

5.2 Constraints on primordial non-Gaussianity from //-distortion fluctuations 17 

6 Conclusions 17 


1 Introduction 

Complete thermodynamic equilibrium between the photons and the electromagnetic plasma exists 
only in the very early Universe in the standard cosmological model. The memory of any event which 
injects energy into the electromagnetic plasma is quickly erased under such circumstances, apart 
from the change in the temperature of the plasma and the photons. This fundamental limit on how 
far and how well we can probe the early Universe using the CMB photons is much further than the 
last scattering surface atz ~ 1100 [1-3] that we observe with the CMB anisotropy experiments such 
as COBE-DMR (Cosmic Background Explorer - Differential Microwave Radiometer) [4], WMAP 
(Wilkinson Microwave Anisotropy Probe) [5], Planck [6], SPT (South Pole Telescope) [7], ACT 
(Atacama Cosmology Telescope) [8] and many others. ' 

An important landmark across which the properties of the Universe change in a significant way 
in a short time is the blackbody surface [9, 10] at z w 2 x 10^. At times earlier than this landmark, 
z > 2 X 10^, the photon creation at low frequencies (x = v/T <s 1, where v is the photon frequency, 
T is the temperature of the electromagnetic plasma, x is the dimensionless frequency and we are 
using natural units with Planck’s constant {h) and Boltzmann constant (^b) equal to unity.) com¬ 
bined with the redistribution of photons over the whole spectrum by Compton scattering maintains 
a Planck spectrum for the photons [9]. The dominant process for the photon creation is the double 
Compton scattering for a low baryon density Universe such as ours [11] with a small contribution 
from bremsstrahlung. At lower redshifts, z < 2 x 10^, the photon creation by double Compton and 
bremsstrahlung becomes inefficient, making the photon number a conserved quantity. Any injection 
of energy is still redistributed by Compton scattering but the inability to create and destroy photons 
means that the equilibrium spectrum is not the Planck spectrum but the Bose-Einstein spectrum (see 
e.g. [12]) with a dimensionless chemical potential (// = -//^/T, where //£• is the thermodynamic 

'See http://http://lambda, gsfc. nasa. gov/links/experimental_sites .cfm for a complete list of CMB ex¬ 
periments 
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chemical potential) that does not vanish in general, 


1 

riBEix) ^ —-■ (1.1) 
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The Bose-Einstein spectrum is a stationary solution of the Kompaneets equation [9, 13]. 

For observational purpose and also for comparison with the y-type distortions [14] it is conve¬ 
nient to define the dimensionless frequency not with respect to the temperature of the Bose-Einstein 
spectrum but with respect to the temperature of a Planck spectrum which has the same number density 
of photons as the Bose-Einstein spectrum, NBE{T,fi) = NpiiTj-ef) ^ {T - ~ 0.456//,Vref = 

v/Tref. In terms of this reference temperature we have 


«BE(.^ref) = 
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where the first term is the Planck spectrum and the second term defines the /i-type distortion in the 
limit of small distortions [9]. Note that at linear order, for small ju, we have // = -heIT ~ -j^ElTrsf- 
We should emphasize that this re-definition is just for convenience of visualization and does not 
affect the analysis of data. This is because the monopole temperature of the CMB is also uncertain, 
at the sensitivity we will be working, and we must fit for the reference temperature simultaneously 
giving us the freedom to choose the reference temperature. We will use the definition from Eq. 1.2 
throughout the rest of the paper with T = 2.725 K as the reference temperature dropping the subscript 

ref. 

The yu-type distortions thus probe the era that lies behind the last scattering surface and are sen¬ 
sitive to any process that injects energy into the electromagnetic plasma [see e.g. 15, 16] between 
redshifts 5 x 10“^ < z < 2 x 10^ [17]. Most processes, owing to the statistical homogeneity of the Uni¬ 
verse, give distortions which are identical in different regions of the Universe. Planck is insensitive to 
the absolute brightness of the sky and therefore cannot detect the constant or invariant component of 
the yu-type distortions. One of the exceptions is the dissipation of primordial perturbations or sound 
waves in the electromagnetic plasma (Silk damping) [18-20] before recombination. The primordial 
perturbations excite standing sound waves in the tightly coupled baryon-photon fluid as they enter 
the horizon [3, 19, 21]. The mean free path of the photons for Thomson scattering on free electrons is 
very small during this time, they can however still diffuse to much larger scales doing a random walk 
in the sea of electrons. The (comoving) diffusion scale (do) and wavenumber (ko) are given during 
fhe radiation dominated era by (including both thermal conductivity and radiative viscosity) [18-20] 
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where R = ^,pb is the baryon energy density, py is the photon energy density, tj is conformal time, 
CTj is the Thomson cross section, ?ie(?7) is the electron number density, //(p) is the Hubble parameter, 
fir is the radiation energy density in units of critical density today, p = 0 is present time and c is the 
speed of light. The diffusion or dissipation scale at z = 5 x 10"* is ~ 46 Mpc“*,/lD ~ 0.14 Mpc 
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and at z = 2 X 10^ it is ~ 1-1 x 10“^ Mpc“^,/}D w 0.55 kpc. The photon diffusion erases the 
perturbations on the diffusion scaies, moving the energy in the perturbations or sound waves into 
the focal CMB monopole spectrum, giving a chemical potential to the spectrum [9, 22-25]. If the 
primordial curvature perturbations are Gaussian, there is same power on the small scales 46 < ko < 
10^ Mpc“^ in different regions of the Universe, apart from the cosmic variance oc 
where < 0.01 Mpc“' is the separation of different regions in the Universe over which we wish to 
study the fluctuations. If the primordial non-Gaussianity of local type is present on these scales [26] , 
the small scale power is correlated with the large scale fluctuations and the ju-distortion has large scale 
anisotropy on the sky with amplitude dependent on the amplitude of the primordial non-Gaussianity. 
The fluctuations in the //-type distortions, and hence the primordial non-Gaussianity on these scales, 
is already accessible with Planck, which has much higher sensitivity compared to the COBE-FIRAS 
for a fluctuating signal [27]. 

We will use the component separation method LIL (Linearized Iterative Least-squares) [28] to 
separate the /i-distortion component from Planck data. The small number of channels and similarity 
of the /i-distortion spectrum with the y-type spectrum means that our map would be dominated by 
contamination from the low redshift y-type distortions and we will only be able to put upper limits 
on the //-distortion fluctuations and primordial non-Gaussianity. In principle, there is additional in¬ 
formation in the intermediate-type distortions [17, 29-31] which is however inaccessible with Planck 
owing to the limited number of frequency channels. The maps and masks, created by one of us (RK), 
are made publicly available at http: //www.mpa-garching .mpg. de/~khatri/muresults/. 


2 /r-type distortions in Planck channels 


We will be using the 4 HFI channels of the Planck (100, 143, 217, 353 GHz). Each of these channels 
has finite bandwidth and we must integrate the intensity spectrum of //-type distortions (2hv^ Ic^n^) 
over each spectral band. In addition it is also convenient to convert the intensity into thermodynamic 
CMB units defined by relating the change in intensity in the spectral band, A/y, to the change in the 
intensity Ip\{T) of the Planck spectrum for a change in temperature ATk^mb’ 


A/y ATKcMB 
“ ^^Kcmb 


dT 

2hv^ xe^ 
{e^ - 1)2 


( 2 . 1 ) 


The units of the maps released by Planck are Kcmb for the 4 HFI channels we will be using. To con¬ 
vert any other spectrum, n{x), to the CMB temperature units we should integrate both the intensity of 
the spectrum of interest and differential Planck spectrum over the transmission profile of the channels 
>v(v) [32] which are available on the Planck legacy archive. 


J w{v)2hv^ I c^n{v)dv 
j w(y)dIpi/dTdv 


( 2 . 2 ) 


For the // distortion n{x) - n^, and similarly for the y-type distortion we substitute n{x) = Uy 
where tiy is the y-type spectrum [14] 
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y to Kcmb 

H to Kcmb 

-5.3337 

-4.1611 

-5.1752 

-2.2528 

-4.7503 

-0.9414 

-4.0309 

-0.2767 

-2.7823 

0.1606 

0.1941 

0.5411 

6.2065 

0.8053 

14.453 

0.9484 

26.332 

1.0413 


Table 1. Conversion factors from y and /r to Kcmb- We have checked that the small difference in the y 
conversion factors with respect to the [32] are due to the change in the transmission profile between the first 
and the second Planck data releases. We have used the later release to calculate the conversion factors. All 
values are for the full channels. 



Figure 1. The spectra of y-type and /r-type distortions in the juKcmb units as seen by Planck channels. We 
will only use the 4 HFI channels, 100 GHz to 353 GHz, in our analysis. 


We give the conversion factors calculated for all Planck channels for both and y-type distor¬ 
tions in Table 1 and the spectra are plotted in Fig. 1 . Our y to Kcmb conversion factors match those 
released by Planck [32] with a small difference because of our using the band transmission profile 
from the second data release. We will exploit the small difference in the two spectra to mask out as 
much of the y-distortion signal from galaxy clusters and groups as possible allowing us to put much 
stronger constraints on the ju-type distortions compared to the constraints we would expect in the 
absence of such a mask. 

3 A yu-type distortion map from Planck HFI data 

The component separation method, LIL, is described in [28] and its application to separate the y- 
type distortion from CMB and foregrounds is described in [33]. The algorithm to create the yu-type 
distortion map is identical to that described in [33] with the y-type distortion spectrum replaced by 
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the //-distortion spectrum given in Table 1. We will only sketch the main features of the algorithm 
below and refer the reader to [28, 33] for details of the implementation. 

LIL is a parameter fitting algorithm and we will fit a parametric model consisting of CMB + dust 
+ //-distortion or CMB + dust -i- y-distortion to the 4 lowest Planck HFI channels at each pixel. The 
spectra of CMB and distortions are of course fixed and the only free parameters are the amplitudes. 
For the dust we fit a gray body dust model with temperature fixed to 18 K but allowing the spectral 
index /3dust to vary between the limited range of 2 < /3dust < 3 away from the galactic plane [see 
28, for exact definitions]. The constraint on the ySdust means that the effective degrees of freedom is 
non-zero even though we fit 4 parameters to 4 data points. This is because away from the galactic 
plane, the main region of interest for us, the signal in 217 GHz channel is too small to constrain jSdust 
and we rarely reach a global minimum inx^ in this parameter direction owing to the hard constraints 
imposed by us. 

An advantage of the parameter fitting approach over the internal linear combination based meth¬ 
ods [34-36] is that we get a quantitative goodness of fit estimate which tells us if the model being fit 
is the correct one. In particular, if we fit a //-type distortion spectrum to the pixels corresponding to 
clusters and groups of galaxies, we will get which would be worse compared to the correct model 
with y-type distortion. We use this difference in x^ between the two models to mask out clusters and 
groups according to the following algorithm 

1. We mask a pixel if -xj, < 

2. We add the requirement that minimum hole size should be 10 pixels at HEALPix [37] resolu¬ 
tion of nside=2048 i.e. if less than 10 contiguous pixels are getting masked then these pixels 
are unmasked. 

3. We also mask all pixels with negative y-type distortion anAx^ >3.8. This is signature of radio 
point sources. 

4. We augment the resulting mask by additional 3 pixels 

This simple algorithm masks most of the clusters and groups that are detected in our y-distortion 
map. In addition we explicitly mask 4.5° radius regions around Virgo and Perseus, 1.5° radius 
regions around all clusters with SIN > 30 in the second Planck cluster catalog [38] and 0.75° 
radius regions around all sources in the second Planck cluster catalog. We finally add the 86% 
CO/point source/galactic mask that was calculated in [33] and is publicly available at http: //www. 
mpa-garching.mpg.de/~khatri/szresults/. Our final minimal mask for the //-distortion anal¬ 
ysis, masking 25.8% of the sky, is shown in Fig. 2. We will augment this minimal mask using 
thresholds on the 545 GHz channel map of Planck to select successively cleaner portions of the sky 
to test the robustness of our results to the dust contamination. 

The probability distribution function (PDF) of the //-distortion map is shown in Fig. 3 for 
different sky fractions /sky- The PDF is dominated by noise and residual y-distortion contamination at 
all sky fractions. For the 62% sky fraction there is also a small amount of contamination from the dust 
(since the masks were designed using the 545 GHz map which is dominated by dust) which becomes 
negligible at /sky = 45% and goes away completely at /sky = 27%. The contribution of the tails to 
the map variance and power spectrum is negligible. This is apparent immediately if we look at the 
same plot drawn on linear scale in Fig. 4. The sharp peak near the small values of // is probably due 
to a small number of clean low noise pixels where the noise dominates over the contamination. The 
skewness towards positive values because of the residual y-distortion contamination is also apparent. 
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74% minimal mask 



Figure 2. Minimal mask for //-distortion analysis masking clusters and groups of galaxies, point sources and 
galactic plane. 


4 Constraints on the //-distortion fluctuations 
4.1 Amplitude of fluctuations 

The simplest constraint we can get on the //-distortion fluctuations is using the total power or variance 
of the map. We should in particular subtract the contribution of the noise. We make half ring //-maps 
from the Planck half-ring channel maps and calculate the half ring half difference (HRHD) //-map 
which gives us the white noise estimate in our full channel //-map. The total map variance cr^ap 
sum of noise contribution, CTnoise and signal //^ms, assuming that the noise is uncorrelated with the 
signal, 


•^map /^rms ^noise’ (4-1) 

where //rms is the root mean squared value of the signal. We are in particular interested in the central 
variance, (//J^g)^ = //^ms - (//)^, where the angular brackets indicate the ensemble average. The 
rms and mean values of // from our map for different sky fractions are given in Table 2. Since our 
algorithm is non-linear, the half-ring half difference (HRHD) map underesfimafes the noise in the full 
channel map. A better, unbiased, estimate of the signal variance would be if we average the product 
of the two half-ring maps. The half-ring maps have uncorrelated white noise which automatically 
cancels in the cross-variance. 


(cr^ap)' ^ 
= 


(4.2) 
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Figure 3. Probability distribution of yu-distortion, P(ij), for different sky fractions, /sky along with the noise for 
/sky = 27% at 10' resolution calculated from the half-ring difference maps. 


where hrl and hr2 refer to the half-ring maps, the ensemble average is over all pixels p, and X refers 
to the cross product of two maps. The central variance is then given by, We 

will use as our upper limit. 

For our smallest mask /sky = 70.7%, the results are clearly dominated by the galactic contami¬ 
nation. The results however vary less rapidly from /sky = 62.3% onward and seem to be converging 
with decreasing sky fractions. Assuming that all the signal is coming from the y-distortion contami¬ 
nation we get an upper limit of < 6.1 x 10“^. Note that since we are averaging over millions 

of pixels, even for the smallest sky fraction, the statistical error is negligible even after taking into 
account the fact that the noise in the neighboring pixels is correlated. 

4.2 Auto and cross power spectrum 

We can also calculate the angular power spectrum of our map, and the cross power spectrum 
with the CMB temperature map, For calculating the power spectra we use publicly available 
PolSpice code [39, 40] which deconvolves the effects of masks [41] in correlation space making it 
very fast. The code also calculates the covariance matrix [42-44] that we use to estimate the errors. 
We use large bins of M = 25 to mitigate to some extent the effects of mode coupling due to our 
very complex mask. We augment our minimal mask using the 545 GHz channel map smoothed to 2 
degree FWHM beam to create the masks with successively cleaner regions of the sky. These masks 
were used in Table 2. For calculating the Cc we apodize the mask with a Gaussian function in pixel 
space by replacing the I 5 in the mask by 1 - exp [-90^/2(0ap)^j for 9 < 0ap = 30', where 6 is the 


- 7 - 
















resolution 10 arcmin 



Figure 4. Same as Fig. 3 but with linear scale on the y axes. 
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18.67 
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1.87 

10.13 

9.95 

7.35 

54.0% 

17.40 

15.15 

2.85 

8.57 

8.08 

6.77 

45.1% 

16.80 

14.77 

3.48 

8.01 

7.21 

6.51 

36.6% 

16.43 

14.48 

3.77 

7.77 

6.79 

6.37 

26.9% 

15.95 

14.14 

3.70 

7.38 

6.39 

6.10 


Table 2. Moments and other properties of the //-distortion map PDF for different sky fractions. The last column 
is the limit on the //-distortion fluctuations calculated from cross-variance of half-ring maps while 
is calculated by subtracting noise estimate from half-ring half difference (HRHD) maps from the full channel 
map. The later slightly overestimates the signal variance. 


distance of the pixel from the edge of the mask. To get an estimate of the power spectrum unbiased 
by noise we will use the cross-power spectrum of the half-ring maps automatically subtracting the 
uncorrelated noise. The auto power spectrum calculated in this way is shown in Fig. 5 for several 
sky fractions corrected for the 10' beam of our maps in addition to the effect of the mask. There is 
a small decrease in the power on small scales when the sky fraction is reduced signifying the level 
of galactic contamination. On large scales, where we have the best constraints, the power spectra for 
different masks are consistent with each other (taking into account the error bars) and we will use 
the /sky - 62.3% power spectrum on large scales to put limits on the //-distortion fluctuations and 
















Figure 5. Auto power spectrum of ju-type distortion fluctuations. The effect of the effective beam of //-map 
(10’ FWHM) and mask has been corrected for. The points for different masks have been offset on the x-axes 
from the center of the Af = 100 bins to make them distinguishable. 


primordial non-Gaussianity. Thus we have the limit on the power spectrum in the i = 1-16 bin, 

i{€ + l)Cf /(27r)|/=2-26 < (2.3 ± 1.0) x IQ-i^ (4.3) 

The cross spectrum of // with primary CMB temperature anisotropies T is given by 

( 4 . 4 ) 

where i is the //-distortion map, either one of the half-ring maps or the full channel map. Similarly j is 
one of the CMB temperature anisotropy maps and is the noise in the respective map. As before we 
can choose the full mission maps and estimate and subtract the noise cross power spectrum estimated 
from the HRHD maps to get = or we can use one of the half-ring //-distortion maps and the other 
half-ring map for T in which case the cross noise term would vanish. We again expect the HRHD 
maps to give slightly biased estimates of the noise because of the non-linearity of our component 
separation algorithm. However on large scales, the cross-spectrum is dominated not by noise but the 
CMB anisotropies and we expect the effect of the noise to be sub-dominant. This can be seen in 
Fig. 6 where we have plotted the cross power spectra between our //-map and the SMICA (Spectral 
Matching Independent Component Analysis) CMB map released by the Planck collaboration [45] 
as well as between our //-map and the CMB map produced by LIL when fitting a y-i-CMB-i-dust 
model [33], so that we have a CMB map where the y-type distortion signal is explicitly removed. All 
power spectra, except one labeled ’hrlx2’, are calculated from the full channel maps and the noise 
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Figure 6. Cross power spectrum of the /U-type distortion fluctuations with CMB temperature anisotropies. The 
CMB temperature anisotropy maps released by Planck collaboration (SMICA) [45] as well as those created by 
LIL in [33] have been used. The effects of the effective beams of the respective maps and masks have been 
corrected for. The points for different masks and map combinations have been offset on the x-axes from the 
center of the bins to make them distinguishable. At high I we use = 100 to make the plot less cluttered. 


power calculated from the HRHD maps has been subtracted. The power spectrum labeled ’hrlx2’ is 
calculated using the first half-ring /t-map and second half-ring SMICA CMB map. The cross-power 
spectrum is clearly dominated by the residual primary CMB anisotropies. This is the same level 
of the signal which is also seen if we cross correlate NILC y-distortion map [46] and the SMICA 
CMB map. Thus it seems that for Planck we are limited not by the sensitivity but by the number of 
channels, the latter limiting our ability to separate the spectral distortions from the primary CMB. 
The SMICA map is much less noisier and less contaminated on large scales because they use more 
channels than the 4 channels used by LIL and we use the cross spectrum with SMICA for 62.3% of 
sky using full channel data to get the limit 


+ l)cf /(27r)|^=2-26 < (2.6 ± 2.6) x lO'^^ K 


(4.5) 


We note that on the large scales all points are consistent while on small scales there is positive 
correlation between the LIL-CMB and maps but negative correlation between the SMICA CMB 
and lu maps. This is just the result of different component separation strategies, in particular the use 
of ^dependent component separation in SMICA. 
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5 Theory of fluctuations of //-type distortions from Silk damping 


The possibility of constraining the primordial non-Gaussianity using the fluctuations of the //-type 
distortions was first presented by [26]. Since then there has been considerable work exploring this 
possibility in detail for many different sources of primordial non-Gaussianity including inflationary 
models with primordial magnetic fields and isocurvature perturbations [47-54]. We will follow a 
new approach to the /i-distortion fluctuations from Silk damping motivated by the interpretation of 
the Silk damping of acoustic modes as mixing of blackbodies [24, 25]. This approach allows for a 
simple but precise treatment of the /i-distortion fluctuations and obviates the need for ’ averaging over 
an oscillation’ and window functions to localize the dissipated energy. We will work in conformal 
Newtonian gauge [e.g. 55] following the conventions for the metric perturbations as defined in [56]. 

In the primordial plasma before recombination, photons and baryons are tightly coupled to¬ 
gether through Thomson scattering. The mean free path of photons is very small but they still diffuse 
to lengths much larger than the mean free path by doing a random walk in a sea of electrons. The 
primordial energy density fluctuations therefore are washed out on this diffusion scale, denoted with 
comoving wavenumber and this process is known as Silk damping [18]. Since the primor¬ 

dial density fluctuations locally had a Planck spectrum because of the thermalization at z > 2 x 10^ 
[9-1 1], the photon diffusion locally just mixes the photons from different blackbodies. The immedi¬ 
ate result of the mixing of blackbodies is y-type distortion [57] which rapidly relaxes to equilibrium 
Bose-Einstein spectrum at z > = 5 x 10^ [9, 10, 15, 17, 58-60]. At smaller redshifts the spectrum 

remains a y-distortion spectrum. The change from y-type distortion to //-type distortion is in reality 
not abrupt but gradual through an epoch of intermediate-type (or /-type) distortions [15, 17, 58, 60] 
where the initial y-type spectrum has relaxed only partially towards the equilibrium Bose-Einstein 
distribution. We have chosen z^ to be near the center of the intermediate-type region and will ignore 
the /-type distortions for simplicity as we do not expect them to affect the constraints from Planck. 
The dissipation of primordial fluctuations not only adds energy but also the photons to the monopole 
component of the radiation. The //-distortion created is then given by the total energy dissipated 
minus the correction from adding photons [9, 10] 


// ^ 



4 AA\ 


(5.1) 


where AE/E is the fractional change in the photon energy density and AN IN is the fractional change 
in the photon number density. This simple equation is just the result of fitting two parameters of 
the Bose-Einstein spectrum, temperature and chemical potential, to two constraints, the given energy 
density and number density of photons. 

Eor Silk damping the above equation implies that [24, 25] 

where x is the comoving position coordinates, 77 is conformal time, and is the CMB anisotropy 

at position x in direction n. We will use bold symbols for vectors with'denoting a unit vector and the 
same symbol in normal font indicating the amplitude of the vector. Note that x and n are independent 
variables (i.e. independent of each other) and describe the 5 of the coordinates in 6 -dimensional 
position-momentum phase space. All terms are functions of rj which we will suppress for simplicity 
unless necessary for clarity. Note that the second term in Eq. 5.1 has resulted only in 1/3 of the 
energy density in the first term going to the //-type distortion [24, 25]. The angular brackets represent 
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an appropriate average that we will describe now precisely. This average, according to our formalism, 
should lead to mixing of blackbodies locally. 

The acoustic dissipation is represented in the Boltzmann hierarchy by the non-vanishing of the 
photon quadrupole, even in the tight coupling limit [see e.g. 56, 61]. This can be understood as fol¬ 
lows. The photons diffusing from different regions/direction to a particular position coordinate have 
different temperatures. The electrons therefore see non-zero anisotropies and Compton scattering 
isotropizes the radiation field seen by fhe elecfron suppressing the anisotropies. Higher the multipole 
order more strongly it is suppressed. The largest anisotropy therefore resides in the quadrupole and 
its dissipation in hydrodynamic terms can be understood as the action of shear viscosity [25]. The 
mixing of blackbodies by Compton scattering is therefore locally just an angular average and this is 
precisely what the angular brackets in Eq. 5.2 represent. Taking the Fourier transform of Eq. 5.2 and 
performing the angular average we get, 


d/i(k) 

d77 




1 

— 0(k',n)0(k-k',n) 
47r 


(5.3) 


where 0(k,n) is the Fourier transform of temperature anisotropy and contains a stochastic part com¬ 
ing from the initial curvature perturbations , !R(k), and a deterministic transfer function which we will 
denote in harmonic space by @i{k). The product of anisotropies in real space gives the convolution in 
Fourier space. We decompose the anisotropies into spherical harmonic coefficients, utilizing the fact 
that at linear order the transfer functions depend only on the scalar product k.n. 

0(k, n) = + l)n(n.k)0f(k)W, (5-4) 

t 


where Vt are the Fegendre polynomials. Substituting in Eq. 5.3 and integrating over n we get 

^ = -2.8^ f ;^y(-l)^(47r)T;jk')kr„Xk^')0r(^')0r(|k-k'|)^R(k')??(k-k'), (5.5) 
drj dT/J (27r)3^ 


where we have used the addition theorem for spherical harmonics, 

m 


(5.6) 


Since we are explicitly doing the angular average of CMB anisotropies locally at position x, our 
energy dissipation is localized by construction and we do not need a window function to localize the 
dissipated energy. 

In the tight coupling limit almost all of the energy is in the monopole and dipole. The tight 
coupling solution for the linear CMB transfer functions is [3, 19, 55, 62-64] 


00 


01 


0, 

A' 


A'(-3 - 0.4/?v) cos(krs)e“^"^^°" 

,( —3 — 0.4/?y) 1 , 2 /r. 2 

A'^-- —- sin{kr,)e-’^ 

V3 

0 if ^ > 2 
0.5 


0.4Ry + l.5 


-0.3 


(5.7) 


where r^{T]) - dr/Cg is the sound horizon, Cg is the sound speed of baryon-photon fluid which af 
z > 5 X 10^ is well approximated by Cg » 1/ V3, kD(r/) is fhe damping wave number [18-20]. These 
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analytic solutions assume that the two metric perturbations (p and ip have the same transfer function 
in the radiation dominated era, 

,Afsin(^'7/V3) - (^ 77 /V3)cos(/:?7/V3) 

= - 

where ' indicates the initial perturbation on super horizon scales in the radiation dominated epoch but 
after neutrino decoupling. These solutions take into account the effect of neutrino anisotropic stress 
in the initial conditions, 0' = -(1 +0.4/?v)<A^ ®() - ~*A72, lA' = KHOARy + 1.5), where Ry - pyjipy + 
Py) w 0.409 is the fraction of neutrino energy density with respect to the total (photons+neutrino) 
energy density in the radiation dominated epoch after the electron-positron annihilation. At the level 
of accuracy needed for us, we can set 3 - 1 - 0.47?v w 3. 

Substituting the tight coupling solutions into Eq. 5.5 and taking the ensemble average of the 
stochastic part, R terms, we recover the results of [24, 25]. It is interesting to note that the sum of 
monopole and dipole terms gives a term oc cos^(krs) - 1 - sin^(krs) = 1 , i.e. the total energy in the 
acoustic oscillations adds up to a non-oscillating constant. This is as expected since the sound waves 
just represent the oscillation of internal energy ( 0 o) and kinetic energy ( 0 i) into each other while the 
total energy of the sound wave must be conserved in the lab frame (at rest with respect to the bulk of 
the fluid). The energy (and therefore the above interpretation) is of course a gauge or reference frame 
dependent quantity and we may interpret the above solutions differently in a different gauge [e.g. 65]. 
In our formalism there is no need to do the average over an oscillation period of the sound wave as 
the total contribution to the //-type distortion does not oscillate. 

Since we are interested in fluctuations of p, we will defer the statistical or ensemble average 
of the stochastic component, R, until the end when we calculate the angular power spectra. At 
high redshifts, z > 2 x 10®, the //-type distortions are suppressed as the photon creation by double 
Compton scattering and bremsstrahlung together with comptonization is able relax the spectrum to 
Planck spectrum, p —> 0. We can take this into account by multiplying the p distortion by the 
blackbody visibility function [9-11], @( 77 ). Integration up to z - z^ » 5 x 10^ then gives us the //-type 
distortions at z^, or the corresponding conformal time 77 ^, 

r''" d r xn / * _ 

//(k, 77^) = - 2.8 d 77 j ^(-1) )k^M(k - k') 

@dk', 77)0r(|k - k'l, T])R(k')R{k - k'). (5.9) 

Precise expressions for 0 are given in [10] but it is well approximated by the simple expression, 
sufficient for our purpose, Q x ^ where ZdC » 2 x 10 ®. 

We will be interested in the squeezed limit, where we want to calculate the fluctuations that are 
much larger than the dissipation scale. The dissipation scale here is given by k' and the fluctuation 
scale by k. We can therefore take the limit k —> 0 in the non-stochastic terms to simplify the above 
expression, 

p{k, 77 ^) - 2.8 £ d 77 §{£ J 0^ £ie + 1 ) 0 ,(k', 77 ) 0 ,(Ik - k'l, r])R{k')R{k - k') 

* -2.53 £ d 77 @( 77 ) A J £^e-^^''^^^"R(k')R(k - k'), (5.10) 

where we have used the tight coupling solution in the last line where the oscillating part evaluates to 
a constant in the limit k —> 0. The time dependent part is now very simple and is trivially integrated 


(5.8) 
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if we approximate @ by a step function that is unity for redshifts in the range 5 x 10"* = < z < Zpi = 

2x10^ and zero otherwise, yielding a final very simple expression for the //-distortion fluctuations 

//(k, 7/^) w 2.53 J ^ - k'). (5.11) 

The fluctuations in // represent fluctuations in the energy density of photons and these energy 
density fluctuations would oscillate in the same way as the primary CMB fluctuations. The // param¬ 
eter that we have defined is however gauge invariant, similar to the dimensionless frequency x - v/T 
and does not undergo oscillations. To be more specific, the adiabatic part of the radiation transfer 
function has no effect on fluctuations of //. The fi fluctuations are therefore frozen in at z - z^- The 
dissipative or Silk damping part of the radiation transfer function would affect the //-fluctuations and 
//-type fluctuations would also dissipate on the diffusion scale which keeps on increasing unfil fhe 
recombinafion. The radiative fransfer of fhe //-fype flucfuafions is fherefore very simple, if is modi¬ 
fied by Silk damping in fhe same way as fhe primary CMB anisofropies unfil recombinafion and affer 
recombinafion fhey jusf free sfream fo us. We can fake info accounf fhe Silk damping very simply fo 
gel fhe //-dislorlion flucfuafions al fhe Iasi scallering surface (LSS), z » 1100, 

//(k, rjLss) = ju(k, (5.12) 

where kD(77LSs) is the damping wavenumber at the last scattering surface which we can evaluate 
numerically for standard recombination history from the well known expression [18-20, 56, 61]. It is 
interesting to note that the dissipation of // fluctuations will result in mixing of Bose-Einstein spectra 
giving a y®^-type distortion different from the y-type distortion coming from the mixing of Planck 
spectra. We can of course trivially calculate this distortion using the same formalism used for mixing 
of blackbodies. For fluctuations around a mean //, we get after mixing a distortion u^be given by 


_ + 1 \ 

2 (^ezc+p _ - 1 j ’ 


(5.13) 


where we have expanded in Taylor series the spectrum - 1). The zeroth order term is 

just the mean Bose-Einstein spectrum, the first order term vanishes since = 0 by definition and 
the second order term gives the lowest order distortion. This distortion would be an effect that is 4th 
order in perturbations, since // is already 2nd order and we will not pursue this interesting possibility 
here. Note that the temperature of the Bose-Einstein spectrum would also vary along with the // and 
will give the usual y-type distortion at lowest order. 

We should note that for modes outside the horizon at the time of recombination, in the Sachs- 
Wolfe (SW) limit [66], the effect of Silk damping is negligible and we can ignore it. 


//^'^(k,77Lss) 


(5.14) 


After the recombination the //-type fluctuations just free stream to us, with a small damping 
at reionization that we will ignore here, and the equation describing it is just the free streaming 
Boltzmann equation. 


rj) ^ - Q (5.15) 

dr] 
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Integrating from t/lss to today t]q we get, 

Mk, 770) = Mk, 

= /^(k, ?7LSs)47r ik{7]o - tjlss)) ^ Jr'm'(n)F;,^,(k), (5.16) 

m' 

where we have used the expansion of the exponential function in spherical harmonics [67]. We can 
now Fourier transform it back to real space and decompose it into spherical harmonic coefficients 
giving, 

/ 13 

(Kno - J]Lss)) F;„,(k)ju(k, ?7 lss) 

/ d^k d^k' 

W - ^Lss)) F;^(k) 

X [g-27c'2/*:D(ZplF _ - k'), (5.17) 

where x is the position of the observer. This is our final expression for /i-distortion anisotropies 
which agrees with the corresponding expression in [26] apart from the trigonometric and the win¬ 
dow functions which are not needed in our formalism. The damping factor for //-type fluctuations, 
g-k^/ko(riLssf ^ can be neglected in the Sachs-Wolfe limit but must be included when using modes in¬ 
side the horizon during recombination {£ > 100). For a statistically homogeneous Universe, such 
as ours is assumed to be, the statistically averaged quantities would not depend on the position of 
the observer. The dependence on x would therefore disappear once we take the statistical average, 
similar to the calculations of the primary CMB anisotropies. To be specific, once we do the statistical 
average of two quantities with wavenumbers ki and k 2 as arguments, such as for the calculation of 
the power spectrum, we will get the factor d^(ki -i- k 2 )e 4 '^i+'^ 2 ) x making the argument of the exponent 
and hence the dependence on the position of the observer vanish [see e.g. 56]. The free streaming 
solution therefore just projects the 3-d fluctuations at the last scattering surface onto a sphere in the 
spherical harmonic basis. This completes our derivation of the evolution of //-distortion fluctuations 
and the formulae are completely analogous to those of primary CMB anisotropies if we make the 
instant recombination approximation and free stream the CMB anisotropies from the LSS to today, 

/ d^k 

(.kim - VLSS)) Y*,^ik)S{k, r]LSs)m), (5.18) 

where S (k, //lss) is the source term that reduces in the Sachs-Wolfe approximation to 

5/7lss) = &o{k, t/lss) + >k(k, t/lss) * (5.19) 

5.1 The auto and cross power spectrum of //-distortion fluctuations 

We will assume the following local model for the primordial curvature perturbations: 

^(x) - ^Rg(x) + ^/nl^g(x)1 (5.20) 

where "Rg is a Gaussian random field and /nl is the non-linear parameter. With this model the cross¬ 
power spectrum between // and CMB temperature oc /^l and the auto power spectrum of //, 
Cf cx 9/25/2^ ^ tnl- 
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We will denote by the angular braekets the statistieal ensemble average from now on. For the 
eross power speetrum we get 

, (5.21) 

where 6 represents the Kroneeker delta funetion and 

10 12 3 r 

= ^3 ^/nl J dkk^ {k{T]o - t/lss)) S {k, T}i^ss)P'RGi^) 

X (■5 22 ) 

where we have used, taking the limit k <^k', 

{K{V!)n{V. - k')^R(k2)) « ^(27r)VNL4(k + k2) [P<RG{k'f + 2P<R^{k')P<R^{k )\. (5.23) 

The first term in the square braekets gives a term of the form J dkk^jt{kx) independent of the CMB 
anisotropies whieh is formally infinite but in the squeezed limit, k ^ 0, should vanish and we diseard 
it ineluding only the seeond term in expression for the The integrals are easy to do for a speetrum 
of form = A7^Q(27r^/k^)(k/ko)”"“^ giving us the final expression. 




.867r/NLA^^(7/o - ?7 lss)' 

Kq 


-/K«s), 


(5.24) 


where. 


ft{n^ tl)^ 

f{(n^ ^ 1) ^ 



^D(Zpl) 
2^(^+l)^^^\kD(z^);’ 


n.- \ 


(^D(Zpl)”‘‘ ^ - kT){Zfjf" 


log 


where T is fhe gamma funetion. 

For the auto power speetrum we get, 
... . ■)2 

cT = 


j dk k^ji (kirjo - ?7lss))^ PnG(k)e 


where we have used 


(5.25) 


(5.26) 


<^R(k'i)^R(ki - k'i)^R(k'2)^R(k2 - k'2)) - + k2) l4P^G(ki)P9iG(kiWk2) + m\,k'^)] , 

(5.27) 


where f(k\ , k'^) represents permutations whieh are independent of ki , k 2 in the squeezed limit, k\,k 2 
0 and would again result in terms in Q of form J dkk^jfikx) whieh should vanish in the limit k ^ 0. 
We have also defined tnl = 9/25/j^j^. The integrals are again analytieally solvable on large seales for 
power law initial eurvature speetrum giving 

10.12^A^ 

=47rTNL —(.m - VLSs)^~""ge(ns), (5.28) 

ko 
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where 


giin^ tl) = 


V^ri 

(¥) 

iri 


4 r| 

[^) 

iri 

(5+2£-n,^ 




gi{n^ = 1) 


1 

ii{e +1) 



^D(Zpl)\ ^ 
kniZfi) / 


(5.29) 


We should point out that we can choose different parameter values for the primordial power spectrum 
on large and small scales in the above expressions. One factor of x the term in square 

brackets in Eq. 5.25 in the expression for corresponds to the small scale perturbations and it is 
possible to choose different A^,kQ, for these factors to take into account, for example, running of 
the spectral index or a break in the power spectrum. The same is true for the term in 

the square brackets in Eq. 5.29 for the 

5.2 Constraints on primordial non-Gaussianity from //-distortion fluctuations 

Using the ACDM cosmology parameters [5, 68], ns = 0.965, = 2.1x10“^, koCz/i) = 46Mpc“\kD(Zp) 

1.1 X 10^ Mpc“',ko = 0.05 Mpc“',77o - t/lss - 13.9 Gpc, we get on large scales 


» 2.4 xIO-'Vnl 

t(.e+ l) ^„ ^ ^ 10 -^Vnl, (5.30) 

2n ‘ 


where we have evaluated the expressions at ^ = 13, the center of our A€ = 25 bin in Planck spectra, 
and the dependence on € on large scales is very mild because is very close to 1. If we choose 
ns - E with the other parameters the same, we get slightly larger signal as expected. 


(ns = 1) = 2.9 X 10-'Vnl K 
(ns = 1) = 2.8 X IO-^Vnl (5.31) 

2n ^ 


Einally we get the constraints on the non-linear parameters using our measurements of the //- 
fluctuations. 


/nl - 
< 

"rNL - 
< 


+ 1 ) 

271 2.4x10-1^ 

10^ 

m+i) cf 

2j: 1.7x10-23 

1.4 X 10 ^^ 


(5.32) 


(5.33) 


6 Conclusions 

We have constructed all sky maps of //-distortion fluctuations from the publicly available Planck HEI 
data. Our maps and masks are made publicly available at http://www.mpa-garching.mpg.de/ 
“khatri/muresults/. Since we have limited number of channels in Planck, especially channels 
not dominated completely by foregrounds, it is not possible to completely separate the //-type and 
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j-type distortions and our /r-type distortions maps are dominated by the y-distortion contamination. 
Nevertheless, it is possible to put interesting upper limits on the anisotropy power spectrum and 
cross spectrum of // with CMB temperature anisotropies. We have presented measured angular auto 
and cross power spectra for { < 1000 which should be taken as the upper limits. 

Our observational constraints on and allows us to constrain new physics which may 
be responsible for spatially varying energy release in the early Universe, 5xl0^<z<2xl0^. 
We apply our results to constrain the primordial non-Gaussianity for extremely squeezed configu¬ 
rations, ks/^L w 5 X 10^ - 10^, where very few constraints exist at present. Our constraints are 
much stronger, for the same scales, than the recent constraints obtained by [53] using only the CMB 
temperature anisotropies. The only other comparable constraints, to our knowledge, come from the 
considerations of formation of primordial black holes from collapse of primordial fluctuations in 
the presence of large non-Gaussianity [69, 70]. Our constraints for the non-linear parameters are 
/nl < 10^,tnl < 1-4 X 10*^ implying that the non-Gaussianity on these extremely small scales, is 
smaller than of order unity. We should note that the mean or average distortion (/r) should be of 
same order of magnitude as the rms, /irms> i-C- a strong non-Gaussianity of local type would also 
increase the average //-distortion compared to the case of no non-Gaussianity, other parameters being 
the same. Thus the constraints from COBE-FIRAS on (//) also constrain primordial non-Gaussianity. 
Our constraints using the Planck data are however much stronger and in particular constrain the non¬ 
linear term /nl^^ to be not larger than the linear term Hq in the primordial perturbation H. With the 
future experiments, such as Pixie (Primordial Inflation Explorer) and EiteBIRD (Eite satellite for the 
studies of B-mode polarization) [71], it will be possible to improve on these limits by many orders 
of magnitude. However, with Planck satellite mission we have already entered a new era of CMB 
spectrum cosmology. 
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